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Abstract 

The persistent current in a mesoscopic ring has a Gaussian distribution with small non-Gaussian corrections. 
Here we report a semiclassical calculation of the leading non-Gaussian correction, which is described by the 
three-point correlation function. The semiclassical approach is applicable to systems in which the electron 
dynamics is ballistic and chaotic, and includes the dependence on the Ehrenfest time. At small but finite 
Ehrenfest times, the non-Gaussian fluctuations are enhanced with respect to the limit of zero Ehrenfest 
time. 


1. Introduction 

The fact that application of a magnetic field in¬ 
duces an equilibrium charge current is at the ba¬ 
sis of the Landau diamagnetic magnetic response 
of metals [lj - For conducting rings threaded by a 
magnetic flux, this orbital magnetic response takes 
the form of a current around the ring, whereas the 
sign of the response may be diamagnetic as well 
as paramagnetic @]. The recognition by Biittiker, 
Imry, and Landauer that this so-called “persistent 
current” continues to exist in the presence of elas¬ 
tic impurity scattering and, hence, should be 
observable in realistic metal samples, initiated a 
surge in theoretical and experimental work on this 
paradigmatic mesoscopic phenomenon in the mid 
1980s and 1990s (3|. Two recent experiments have 
revived the interest in persistent currents SS0- 
The magnitude of the measured mean square cur¬ 
rent is in excellent agreement with the original the¬ 
oretical predictions for disordered metal rings ;9]. 
Earlier experiments had confirmed the existence of 
the persistent currents Ed, El, but a quantitative 
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verification of the theoretical estimates was not pos¬ 
sible. 

Whereas disorder is unavoidable in metal rings, 
persistent currents were also investigated in semi¬ 
conductor heterostructures, for which the electron 
motion is ballistic 12], The most pronounced differ¬ 
ence between ballistic and disordered-diffusive rings 
is the possible existence of short periodic electron 
trajectories in the former, for which the persistent 
current essentially follows the behavior of ideal one¬ 
dimensional rings without potential scattering 13]. 
Such short trajectories may dominate the magnetic 
response, even if the classical dynamics in the bal¬ 


listic conductor is chaotic [1J, 15, 16, 17]. 


An interesting case arises if the ballistic conduc¬ 
tor has a chaotic classical dynamics, but without 
short periodic trajectories encircling the magnetic 
flux [181 ]. Examples of such a situation are, e.g., 
a ballistic ring with disc-like scatterers, referred to 
as a “Lorentz gas”, or a collection of chaotic cavi¬ 
ties arranged in a ring. Without short periodic tra¬ 
jectories, differences between the ballistic chaotic 
conductor and its disordered counterpart are much 
more subtle, related to the “Ehrenfest time” r E [l9j . 


r E = ^ In kL, 


( 1 ) 
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where A is the Lyapunov exponent of the classical 
dynamics, k is the wavenumber, and L a charac¬ 
teristic classical length scale. Being the time re¬ 
quired for two classical trajectories a quantum sep¬ 
aration 1/k apart to acquire a classical separation L 
under the influence of the chaotic classical dynam¬ 
ics, te characterizes the threshold between classical- 
deterministic and quantum-stochastic dynamics in 
ballistic structures. Ehrenfest-time-related effects 
have been considered for equilibrium properties of 


chaotic quantum dots 00 


turn transport in o pen systems 



, and for quan- 

M 00(27, 


0 0 0 dll . 0. Id .'ll. 0, but not for persistent 
currents in a ring geometry. 


In the present article we report a study of the 
Ehrenfest-time dependence of the mesoscopic fluc¬ 
tuations of the persistent current in ballistic rings 
in which the classical electron motion is chaotic 
and, after appropriate coarse graining, diffusive. 
We consider a grand canonical ensemble, and as¬ 
sume that time-reversal symmetry in the ring is 
broken by an applied magnetic field. In a ballis¬ 
tic ring, mesoscopic fluctuations of the persistent 
current are induced by variations of the chemical 
potential /u; no disorder average is taken. Differ¬ 
ences between ballistic-chaotic conductors and their 
disordered counterparts appear through a depen¬ 
dence on the Ehrenfest time 7 e for the ballistic- 
chaotic case, whereas te plays no role in the case of 
a disordered conductor. As we show below, no te- 
dependence is found on the level of the two-point 
correlation function (/(</>i)/(</> 2 )} of the current dis¬ 
tribution; Only the connected three-point correla¬ 
tion function K(<j)= (/(<(>i)/(</> 2 )/(</> 3 )}c, 
which describes deviations from the Gaussian dis¬ 
tribution, shows a dependence on the Ehrenfest 
time in the case of a ballistic conductor. (Here (j) 
is the flux threading the ring, in units of the flux 
quantum /ic/e; The subscript ‘c’ refers to the ‘con¬ 
nected average’, ( abc) c = (abc ) — ( ab)(c} — ( bc)(a ) — 
(ca){b) + 2(a)(6)(c).) 


Below, in Sec. [2] we describe the starting point 
of our theoretical approach, Gutzwiller’s trace for¬ 
mula, and the semiclassical approximation. A cal¬ 
culation of the two-point correlation function is pre¬ 
sented in Sec. [3j and the three-point correlator is 
discussed in Secs. [4] and [5] We conclude in Sec.[G] 


2. Persistent current from Gutzwiller’s trace 
formula 


Starting point of our calculation of the persistent 
current I is the thermodynamic relation 


e dfl 
h d(j) ’ 


( 2 ) 


where the thermodynamic potential at temperature 
T and chemical potential /z, 


n = -T J de ln(l + e- (e -^ /7 >(£), (3) 


is expressed as an integral of the density of states 
v(e). Following previous works on persistent cur¬ 
rents in ballistic chaotic conductors Ibllfbl [TtH. 17 1, 
we use the Gutzwiller trace formula |35| to express 
the fluctuating contribution to the density of states 
as a sum over periodic orbits a on the energy shell 


i/(e) = eY J A a t°y s ^ /h . (4) 

a 

In this expression, the label a represents a periodic 
orbit with primitive period and period t a = , 

where m is the repetition number. Further S a (e) 
is the classical action of the orbit a and A a the 
stability amplitude of the orbit, 

A a = [det((M°) m - 1)]" 1/2 (5) 

where is the stability matrix of the primitive 
orbit a 0. 

We now specialize to a two-dimensional system 
threaded by a flux <f> = </>/ic/e. Considering energies 
e near the chemical potential /z, the action S a (s, (f>) 
can be written 


S a (e, (/)) = S a (n, 0) + 27 T(j)hn a + (e — (6) 

where n a is the winding number of the trajectory a. 
Below we will write S a as short-hand notation for 
S a (iJ,, 0). Substituting the Gutzwiller trace formula 
for the density of states iq taking the derivative to 
cj), and performing the integration over £, one finds 


I = - 


xe 


E 




2nh ^ t a sinh(7rf a T /H) 


( 7 ) 
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Upon separating the current into Fourier compo- phase space point at time t , while having passed n 
nents, times around the flux, we conclude that 


I = J2 Ine 2 ™ 4 ’, 


with I n = I*_ n , one then arrives at the result 


len 

ln = 


7T Tt° 


( 8 ) 


( 9 ) 


2irh ^ t a sinh(7 xt a T/K) 
x (A a e lS °/ h 6 na}n + A* a e- iS “/ h 5 na ,- n ) . 


3. Mean square current 

We now calculate the mean square (/„/_„} for 
the case that time-reversal symmetry in the ring is 
broken by an applied magnetic field. The leading 
contribution to ( I n I- n ) comes from diagonal con¬ 
tributions, 


(Inl-n) = 


2^2 


2 e n 

(2irh) 2 


( 10 ) 


v (tt T)\ty ta y 4 

a sinh 2 (7 Tt a T/h) 


The factor two in the numerator comes from the two 
terms in Eq. m, which give equal contributions to 

(irJ-n). 

In order to perform the trajectory sum in Eq. 
CD, we use a method proposed by Argaman, Irnry, 
and Smilansky 37]. The summation over classical 
trajectories is expressed as an integral over the en¬ 
ergy shell Q. Introducing a phase space coordinate 
/x, and denoting with fi(t) the phase space coordi¬ 
nate obtained by following the classical time evolu¬ 
tion for a time t, starting at /x, one has 


]t° a \A a \ 2 S ncttn 6{t - t a ) 


( 11 ) 


'Q 




{Inl-n) = 


dt ■ 


(nTf 


2ir 2 h 2 J t smh 2 (ntT/h) 

d/xp(/x,/x,f|n). (13) 


Here we neglected the contribution from orbit rep¬ 
etitions, which is a standard approximation in this 
field, since the non-primitive orbits at a given pe¬ 
riod are exponentially outnumbered by primitive 
orbits with the same period. 

For a two-dimensional ring of circumference L 
with diffusive electron dynamics, one has 


P(V, V,t\n) = 


L e ~(nL) 2 /4Dt 

Q \/ AirDt 


(14) 


where Q = 27t/ith is the volume of the energy shell, 
th being the Heisenberg time, and D the classical 
diffusion constant. One then arrives at the result 


{Inl-n) = 


2 2 
e n 



where 


TL = -D 


(15) 


(16) 


is the time required to diffuse arouud the ring. This 
is the same result as what one obtains for a disor¬ 
dered metal ring [§;,;§]. In the limit of zero temper¬ 
ature, Eq. CD simplifies to 


(Inl-n) = 


6e 2 


TT 2 n 3 Tj 


(17) 


For high temperatures, T K/tl , the integration 
can be performed using the saddle-point method 
and gives 


where n(/x, t) is the number of times the trajectory 
starting at the phase space point fi winds around 
the flux in the time t. The factor t° a arises, be¬ 
cause each trajectory is weighted by a factor t° a 
upon performing the phase space integration |36I |. 
Upon identifying 

~ M)<5n( M ,t),n = p(v, /x, t\n), (12) 

as the classical probability density that a particle 
starting at phase space point /x is found at the same 


(Inl- n ) = |n|T 2 e-Hv^^A, (18) 

up to corrections that are small in the limit T 
H/t l . 

The main result of this section is that the two- 
point correlation function is the same for a ballistic 
chaotic ring and for a disordered metal ring, pro¬ 
vided the coarse-grained classical dynamics in the 
ring is diffusive. The Ehrenfest time te of Eq. ID 
has not entered into our considerations. 
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Figure 1: Trajectories a (dashed), o! (solid), and cx" (solid) 
contributing to the expectation value The 

trajectory a has a “figure eight” structure and is piecewise 
equal to ex' or a/ 7 , up to quantum uncertainties. The differ¬ 
ences between a on the one hand and ex' and ex" on the other 
hand are maximal in the TE-long “encounter” (dark shaded 
part) in which cx “switches partners”. The points /i and y! 
denote the points in phase space where the self-encounter 
starts and ends respectively, (a) The duration of the self¬ 
encounter, te, can be shorter than both t a / and t a n. (b) 
The self-encounter can wrap w times around the shorter of 
the orbits ex' and ex". Here cx" is assumed to be the shorter 
one, and we illustrated the case ofw = 1. 


4. Third cumulant: Off-diagonal contribu¬ 
tion 


The second moment of the current distribution 
could be calculated by considering diagonal terms 
in the trajectory sum only. For the calculation 
of the connected expectation value 
the Fourier transform of which gives the connected 
three-point function K((f>\, </> 2 , <(> 3 ), one needs to go 
beyond the diagonal approximation. It is at this 
point, that the Elirenfest time te enters into the 
calculation flflj . 

The semiclassical calculation of a connected 
three-point function of the density of states (which 
contains all essential information for the three- 
point function of the persistent current) was per¬ 
formed by Heusler and coworkers for the case of 


a chaotic quantum dot [38j, [39|, building on previ¬ 


ous developments of the trajectory-based semiclas¬ 


sical formalism by Sieber and Richter [ioj]. Fol¬ 
lowing Refs. [38|,[39|, the dominant contribution to 
{InIm,I-n-m)c is given by a summation over “tra¬ 
jectory triplets”. These trajectory triplets consist 
of a trajectory a which contains a small-angle self¬ 
encounter, so that it effectively has a “figure-eight” 
structure, see the dashed path in Fig. [T}r. The other 
two trajectories a' and a" (solid paths) are differ¬ 
ent and piecewise equal to one of the loops of a, up 
to a quantum uncertainty. Once the trajectory a 
is specified, the other two trajectories are uniquely 
determined, so that a' and a" need not be summed 
over separately. The periods of the three trajecto¬ 
ries in Fig. Q] are related as 

toe — to,' — t toe" (19) 

and for the stability amplitudes one finds 


A n = An'An 


( 20 ) 


Accounting for the various ways in which the trajec¬ 
tories can be combined, one obtains the expression 


{Inlml—n—m) c 

2 ie 3 nm(n + m) 


^ ' eir 


tT\A d 


(27 t) 3 H 3 X—j sinh(nt a T/h) 

(tt r) 2 Ree i ( s «- s »'- s «'")/ fi 
sinh(7rt a /T/?i) sinh(7 xt a "T/h) 


X ^ n a r ■ rl < ^n, y " ,m “b ,n&', 

“b 


n n // ,—m—n 


J n n // ,—m—n) 


( 21 ) 

The summation over trajectories is now per¬ 
formed using the method of Refs. SHI and their 
extension to the systems with diffusive classical dy¬ 
namics [ 32 }. The oscillating factor in the numer¬ 
ator of Eq. m suppresses contributions from all 
trajectory triplets for which the action difference 
AS = S a — S a ' — S a " between a on the one hand 
and a! and a" on the other hand is not at most of 
order h. Since the action difference is related to the 


duration r enc of the small-angle encounter 41, l42j . 
|A<S/ft| ~ kLe ~ XTetl °, one finds that only trajectory 
triplets with r enc = te contribute to the summa¬ 
tion. Proceeding as in Ref. 32] one then finds that 


^2 \A a \ 2 cos(AS/h)S naltn S na „ jm 5(t a ' -ti) ( 22 ) 


5{t a " — t2 ) = ^2 


w>0 


1 d 

r H <9 t e 


A(ti,t 2 ;T E ; w), 


where the function F(ti,t- 2 ',TE',w) is the properly 
normalized probability for the trajectory configu¬ 
ration to occur. It depends on the number of times 
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w the encounter winds around the shorter of the 
orbits a' and a", 

F(ti,t2\T E -,w) (23) 

= Q J dpdp'p(p, p!,ti — r E \n — wm) 

x p(p',p,T E - wt 2 ) 
x p(p, p!, t 2 — r E + wt 2 \m) 

if ti > t e and t e /(uj + 1) <t 2 < t e /w, i.e., if t 2 is 
the duration of the shorter orbit and the encounter 
wraps w times around it, 

F(h,t 2 ;T E ] w) (24) 

= Q J dpdp'p(p, p! ,t 2 — TE\m — urn) 

X p(p',p,T E - wti) 
x p(p, p',h - r E + whin) 

if t 2 > r E and r E /(w + 1) < h < t e /w, i.e., if h is 
the duration of the shorter orbit and the encounter 
wraps w times around it, and 


F(h,t 2 -,T E ]w) = 0 


(25) 


if both h and t 2 are smaller than r E , in which case 
the figure-eight configuration of Fig. |T] is not possi¬ 
ble because the trajectories a ' and a" are identical 
and a is a non-primitive orbit. Figure |T}r shows an 
example of an orbit configuration with w = 0; An 
example with w = 1 is shown schematically in Fig. 

ID). 

Denoting the distance around the ring’s circum¬ 
ference between the phase points p and p! point by 
x, we have 

T p—(nL—x) 2 /4Dt 

p{p,p',t\n) = 


p(p',p,t) = 


Q y/AirDt 

L e -x 2 /ADt 


(26) 


Q \JAirDt 

where we do not impose a bound on x to account 
for the possibility that the encounter itself winds 
around the ring. Substituting these explicit expres¬ 
sions for the probability densities, one finds 

tl 


F(h,t 2 ;T E ;w) = 


ATTy/a{t 2 ,h]T E ;w) 

_ t l (n 2 t 2 -2mnT E 4-m 2 (t 1 -m(u) + l)t 2 +2wT E )) 

X e 4o-(t 2 1*i; te; w ) 

if ti > r E and t e /(w + 1) < t 2 < t e /w, 


(27) 


F(h,t 2 ;TE;w) = 


tl 


4tt cr(ti, t 2 ; t e ; w) 


(28) 


Tt fn. i ll-2».T E +. 1 (l 2 -»(»+l)l 1 +2»r i ;)) 

4er(t 1 


if t 2 > r E and t e /(w + 1) < h < t e /w, and 

F(h,t 2 ;T E ;w) = 0 (29) 

if ti < r E and t 2 < r E . Here 

<r(ti,t 2 ‘, r E ; w) = hh - r| - w(wh + h - 2r E )ti. 

We note that F is continuous at fi j2 = t e /w with 

w = 1 , 2 ,.... 

We first perform the remaining integration over 
ti and t 2 in the limit te C tj.. In this limit, it is 
sufficient to consider the case w = 0 only, and we 
may take the limit t e /tl —> 0 after differentiation 
to r E . We find 

^ ^ |A a | C0s(A<S/tl)(5 na , ,ndn a /f ,md(t a ' tl) 

ot 

X 6(tail - t 2 ) 

_. mnr L -m 2 T L /4t 2 -n 2 T L /4 : t 1 /onx 

“ Sttth^) 3 / 2 • ^ 

Performing the remaining integrations over ti and 
t 2 in the limit of zero temperature then gives 


|A a | 2 cos(A<S) 

/ J ~ ~ ~ On a /,nO: 


t at a't cl" 

— f 

9 Jn.mi 


n n ii ,m 


(31) 


with 


fn,m = sign (mn) 


\m\ 2 + 4|m||n| + \n\ 2 
to 2 7i 2 (|to| + |n|) 4 


Hence, in the limit r E -C tl and at zero tempera¬ 
ture, one finds 

(d-ndml—n—rn) 

— (I—nd—mln+m) 

3 e 3 mn(n + to) 


i7r 3 r|r H 


(32) 


X (fn,m 4“ fn, — m—n 4“ fm, — m—n)• 


This result is the same as that was found previously 
for disordered metal rings |43]. 

Including a finite Ehrenfest time r E into the zero- 
temperature calculation leads to a modification of 
the coefficients f n ,rm which now acquire a depen¬ 
dence on t e /t e . We were not able to perform the 
integrations over t\ and t 2 in closed form at finite 
Ehrenfest time, but the integrals can be evaluated 
numerically. The Ehrenfest-time dependence of f\,\ 
and f\- 2 is shown in Fig. Ot and the resulting cu- 
mulant (I±IiI- 2 ) is plotted in Fig. [2}r in units of 


x e 
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Figure 2: (a) Ehrenfest-time dependence of the dimension¬ 
less coefficients /i,i (solid) and /i ,—2 = f— 2,1 (dashed), 
(b) Left: The magnitude of the third cumulant ( I 1 I 1 I— 2 ) 
in units of e 3 / in 3 shown for zero temperature (solid 

line) and Ol = 10 (dashed line). Right: The same zero- 
temperature result shown for small te /tl (solid line) and the 
small -te expansion of Eq. (133K (dotted line), (c) Tempera¬ 
ture dependence of the third cumulant at t-^/tl = 42/999 ~ 
0.042, which is close to its maximum at Ol = 0. 


e 3 /i7r 3 r|TE (solid line). Two remarkable observa¬ 
tions are in place: (i) For moderate but still small 
values of te/tl, the inclusion of a finite Ehrenfest 
time leads to a rather significant enhancement of 
the non-Gaussian fluctuations, (ii) For larger val¬ 
ues of te/tl the cumulant changes sign. 

In the physically relevant limit of small te/tl, 
contributions with t\ < te or t -2 < te are expo¬ 
nentially small in the large parameter tx/te, so 
that it is sufficient to consider the integral for times 
ti ,2 > te, for which one can take Eq. (1271) (1751) with 
w — 0. The result of a series expansion in the small 
parameter te/tl then yields 

f F fc (n,m)(r E /r L ) fc 

n,m (rnn) k+1 |mn|(|m| + |n|) 2fc + 4 ’ 


where the first three coefficients Fk(n,m) are 

Fo(n,m) = m 2 +4\mn\ + rr, 

Fi(n,m) = 16(m 4 + 6|m 3 n| 

+ 15m 2 n 2 + 6|mn 3 | + n 4 ), 
^(n, m) = 240(m 6 + 8|m 5 n| 

+ 28 m A n 2 + 56|m 3 n 3 | 

+ 2877T 2 n 4 + 8|mn 5 | + n 6 ). 

In Fig. [2}o (right plot) we show the cumulant 
(I 1 I 1 I- 2 ) resulting from this second-order expan¬ 
sion (dotted line) together with the full numerical 
solution (solid line). 

For temperature T> H/tl we can perform the 
integrals over t\ and £2 using a saddle-point ap¬ 
proximation. In the limit of small Ehrenfest times 
te <C tl one finds 


(7n7 m /__ n—rn) 

e 3 T 3 TLmn(n + m) 

IT^H 3 

X (*?m,n T 9n,—m—n T 9m,-m-n)i 


(34) 


with 6l = 2irTTL/h 1 and 


g m , n = sign (?rm)e- (|m|+|n|) ^ r . (35) 


We can extend this result to finite Ehrenfest time, 
yielding 


3m,n(r E ) = sign(m?z)e- (|m|+l " l)VFr e 2sign(mn)eE , 

(36) 

with 0e = 2irTTE/h -C 0l- Again, as in the zero- 
temperature case, at finite temperatures a finite 
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Figure 3: Trajectories a , a', and a" contributing to the 
expectation value ( I\IiI— 2 ). The trajectory a is a non¬ 
primitive orbit, consisting of a twofold repetition of a? = a". 
The trajectories a' and a" each wind once around the ring. 


Ehrenfest time can actually lead to an increase of 
the third cumulant. 

For general temperature, we have to evaluate 
Eq. (Eli) numerically. In the left panel of Fig. [SJd 
we show the cumulant (I 1 I 1 I- 2 ) as a function of 
te/tl at 9l = 10 (dashed line). The qualita¬ 
tive behavior of the cumulant is the same as at 
zero temperature: Initially, its magnitude increases, 
until it reaches a maximum at small but finite 
te/tl ■ At longer te it decreases again, eventu¬ 
ally changing sign at te/tl ss 0.2. I 11 Fig. [3^ we 
show the temperature dependence of (I 1 I 1 I- 2 } at 
te/tl = 42/999 ~ 0.042, which is close to the po¬ 
sition of the zero-temperature maximum. 

5. Diagonal contribution to the third cumu¬ 
lant 

In addition to the off-diagonal contributions to 
the third cumulant that were discussed in Sec. 01 
there are also diagonal contributions that involve 
orbit repetitions. Such contributions are usually ne¬ 
glected in a semiclassical analysis, because they are 
suppressed with a factor e^ At , where t is the period 
of the (primitive) orbit and A the Lyapunov expo¬ 
nent. Since the period of typical orbits that encircle 
the ring tl 1/A, one argues that such contribu¬ 
tions can safely be neglected. However, diagonal 
contributions do not involve the inverse phase-space 
volume, so that they lack the factor tl/th that sets 
the smallness of the off-diagonal contributions such 
as Eq. (l32l) . 

The leading such diagonal contribution (for given 
to and n) requires the two short orbits a' and a" 
to be ?r-fold and m-fold repetitions of a primitive 
orbit «o with unit winding number, whereas the a 
is the of the (to + n)-fold repetition of the same 


orbit, with to, n > 0. The case n = m = 1, for 
which a' = a" = a 0 , is illustrated in Fig. [3] For 
such a diagonal contribution one has a single sum 
over orbits, 


(l 


IT \ d 

n ± m ± —m—n/c 


(37) 


47t 3 7i 3 


J2 C ^m(t ao )S nao ,l\A ao \ 2n+2rn , 


where we used that we may set A a i = A” o , A a n = 
A™ 0 , and A a = A™ +n . We abbreviated 




(ttT) 2 


smh(nirtT / h) sinh(m7rfT/ K) 
7 rT 

y/ _ 

sinh[(m + ri)wtT/h ] 


(38) 


For uniformly hyperbolic dynamics one has 
|A ao | 2 = l/2sinh(Af ao ) ~ °~ xt 
This then gives 


e~ xta o for Xt ao > 1 441. 


(J-n Im I— m—n ) c 


ie w 


47T 3 /l 3 


^ ^ Cn,ra(^a 


x^ ao ,i|A ao | 2 e-(™+”- 1 ) At “o. 


(39) 


Performing the remainder of the calculation as in 
Sec. [3] one finds that the diagonal contribution to 
the third cumulant of the persistent current reads 


(Iniml— n—rn) c 
„-„3 


47T 3 7i. 3 




I TL r ~T L /it 

47 Tt 


x e 


— (m-\-n—l)Xt 


(40) 


For temperatures T <C one then finds an 

essentially temperature-independent diagonal con¬ 
tribution to the third cumulant of the persistent- 
current fluctuations, 


IT T T \ d 
\ J n ± m- L —n—m / c 


(41) 


e 3 A 3 / 2 (to. + n - I ) 3 / 2 


7 t 3 t 3 / 2 rnn(m + n) 


At zero temperature all diagonal contributions 
from orbit repetitions are smaller than the off- 
diagonal contributions if the condition 



e yJ (m+n- 1)At£ 

(At l ) 3 / 2 


(42) 


is met. This condition can also be rephrased 
in terms of the Ehrenfest time te, using te ~ 
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A 1 1ii(th/tl), as 


t e < 


(n + m — 1 )t l 
In (t h /tl) 


(43) 


Given the intrinsic smallness of the Ehrenfest time, 
this condition is easily met. However, to see a non¬ 
trivial Ehrenfest-time dependence of the persistent 
current fluctuations, the Ehrenfest time needs to be 
a finite fraction of tl, see Sec. [H and the condition 
( 1551 ) effectively limits the applicability of the results 
as shown in Fig. [2] to the range te/tl <C 1, where 
the expansion ( 1551 ) is valid. 

At finite temperatures T > H/tl the off-diagonal 
contribution is strongly suppressed and the diago¬ 
nal contribution quickly takes over. This reflects 
the large difference in typical orbit durations for 
the off-diagonal and diagonal contributions: For the 
off-diagonal contribution, the typical orbit dura¬ 
tion ~ tl at zero temperature, so that temperature 
starts to suppress this contribution for T > K/tl- 
The typical duration of orbits contribution to the 
diagonal contribution is ~ y/ix/A, which explains 
the relative insensitivity of this contribution to tem¬ 
perature. 


6. Discussion and conclusion 


The distribution of the persistent current in 
a mesoscopic ring is Gaussian, with small non- 
Gaussian corrections. Here we have presented 
a semiclassical calculation of the leading non- 
Gaussian correction, described by the three-point 
correlation function K = (/(</>i)/(/> 2 )d(/> 3 ))c- In 
agreement with previous work for disordered metal 
rings 43, 4fj, we found here that at small temper¬ 
atures K ~ e 3 /r£rH = e 3 /gr 3 , where g = te/tl is 
the dimensionless conductance of the ring, tl the 
diffusion time, and th the Heisenberg time. The 
semiclassical approach also contains information on 
the role of the Ehrenfest time in such a ring, and 
we showed that for small but finite te/tl the mag¬ 
nitude of the non-Gaussian corrections is enhanced 
by a numerical factor, before it is suppressed in the 
limit of large te/tl . 

The fact that the three-point correlation function 
initially increases with increasing Ehrenfest time is 
remarkable, since a finite Ehrenfest time usually 
suppresses quantum interference effects. However, 
it is not without precedent: The conductance fluc¬ 
tuations in a chaotic cavity are Ehrenfest-time in¬ 
dependent [271, [26(, whereas the conductance fluc¬ 


tuations in a quasi-one dimensional Lorentz gas are 


larger in the limit of large Ehrenfest time than in 
the limit of zero Ehrenfest time [32|. The same ap¬ 
plies to the variance of the current pumped through 
a chaotic cavity with a periodic modulation of its 
shape jjifij. Just as in the present case, the conduc¬ 
tance fluctuations or the variance of the pumped 
current contain contributions from closed loops, 
and it is this type of correction that can in prin¬ 
ciple be enhanced by Ehrenfest-time corrections. 

We have also identified a second semiclassical 
contribution to K, which involves a diagonal or¬ 
bit sum with non-primitive periodic orbits. Non¬ 
primitive orbits are usually neglected in semiclassi¬ 
cal approaches, because their contribution is expo¬ 
nentially suppressed in comparison to contributions 
from primitive orbits. Nevertheless, for the three- 
point correlation of the persistent current, such di¬ 
agonal contributions become dominant in the limit 
of large Ehrenfest times and/or temperatures. A re¬ 
lated (but not identical) effect appears for conduc¬ 
tance fluctuations in a chaotic quantum dot, where 
“classical fluctuations” become important for large 
Ehrenfest times 


27]. 


The small magnitude of the non-Gaussian fluc¬ 
tuations turns its measurement into a considerable 
challenge, even with state-of-the-art techniques |hj, 
15]. For disordered metal rings, the conditions for 
measuring non-Gaussian corrections to the distri¬ 
bution are most favorable if its dimensionless con¬ 
ductance g is not too large, since one needs to av¬ 
erage over at least ~ g 2 statistically independent 
samples to be able to distinguish the three-point 
function from the Gaussian (second-order) fluctua¬ 
tions For Ehrenfest-time-related corrections to 
become relevant merely making g small is not the 
solution, since a small dimensionless conductance g 
also implies a small te. 
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